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Modifications of general relativity provide an alternative explanation to dark energy for the ob- 
served acceleration of the universe. We calculate quasilinear effects in the growth of structure in 
f{R) models of gravity using perturbation theory. We find significant deviations in the bispectrum 
that depend on cosmic time, length scale and triangle shape. However the deviations in the reduced 
bispectrum Q for f{R) models are at the percent level, much smaller than the deviations in the bis- 
pectrum itself. This implies that three-point correlations can be predicted to a good approximation 
simply by using the modified linear growth factor in the standard gravity formalism. Our results 
suggest that gravitational clustering in the weakly nonlinear regime is not fundamentally altered, at 
least for a class of gravity theories that are well described in the Newtonian regime by the parameters 
Gcfi and This approximate universality was also seen in the N-body simulation measurements 

of the power spectrum by Stabenau & Jain (2006), and in other recent studies based on simulations. 
Thus predictions for such modified gravity models in the regime relevant to large-scale structure 
observations may be less daunting than expected on first principles. We discuss the many caveats 
that apply to such predictions. 



PACS numbers: 98.65.Dx,95.36.+x,04.50.+h 



I. INTRODUCTION 

The energy contents of the universe pose an interesting 
puzzle, in that general relativity (GR) plus the Standard 
Model of particle physics can only account for about 4% 
of the energy density inferred from observations. By in- 
troducing dark matter and dark energy, which account 
for the remaining 96% of the total energy budget of the 
universe, cosmologists have been able to account for a 
wide range of observations, from the overall expansion of 
the universe to the large scale structure of the early and 
late universe [1]. 

The dark matter/dark energy scenario assumes the va- 
lidity of GR at galactic and cosmological scales and intro- 
duces exotic components of matter and energy to account 
for observations. Since GR has not been tested indepen- 
dently on these scales, a natural alternative is that GR 
itself needs to be modified on large scales. This possibil- 
ity, that modifications in GR on galactic and cosmolog- 
ical scales can replace dark matter and/or dark energy, 
has become an area of active research in recent years. 

Attempts have been made to modify GR with a fo- 
cus on galactic [2] or cosmological scales [3-5] . Modified 
Newtonian Dynamics (MOND) and its relativistic ver- 
sion (Tensor- Vector-Scalar, TeVeS) [2] attempt to ex- 
plain observed galaxy rotation curves without dark mat- 
ter (but have problems on larger scales). The DGP 
model [4], in which gravity lives in a 5D brane world, 
naturally leads to late time acceleration of the universe. 

Adding a correction term f{R) to the Einstein-Hilbert 
action [3] also allows late time acceleration of the universe 



to be realized. 

In this paper we will focus on modified gravity (MG) 
theories that are designed as an alternative to dark en- 
ergy to produce the present day acceleration of the uni- 
verse. In these models, such as DGP and f{R) models, 
gravity at late cosmic times and on large-scales departs 
from the predictions of GR. By design, successful MG 
models are difficult to distinguishable from viable DE 
models against observations of the expansion history of 
the universe. However, in general they predict a differ- 
ent growth of perturbations which can be tested using 
observations of large-scale structure (LSS) [6-20] . 

In this paper we consider the quasilinear regime of 
clustering in which perturbation theory calculations are 
valid. We explore what f{R) modifications to gravity 
predict about the behavior of the three-point correlation 
function. In §11 we outline the particular type of f{R) 
gravity model we will be using for our calculations. In 
§111 we introduce the fundamentals of perturbation the- 
ory and in particular how it applies to modified gravity. 
In §IV we focus on second order corrections and in par- 
ticular the Bispectrum. In §V we present our results and 
compare with other studies of the nonlinear regime. In 
§VI we discuss the implications for observations. 



II. THE f{R) MODIFIED GRAVITY MODEL 

In general f{R) models are a modification of the 
Einstein-Hilbert action of the form: 
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R + f{R) 
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where R is the curvature, = SttG, and Cm is the mat- 
ter Lagrangian. A major issue with gravity modifications 
has been that while they are successful at explaining the 
acceleration of the universe, they also tend to fail to com- 
ply with Solar system (very small scales) observations. 
Recently, though, the so called Chameleon mechanism 
was found [21] that makes it possible to overcome this 
problem. Significant attempts have been made to in- 
clude Chameleon behavior [22] in DGP theories as well. 
One example of an f{R) model that exhibits Chameleon 
behavior was constructed by Hu and Sawicki [23]. The 
functional form of /(i?) there is derived from a list of 
observational requirements: it should mimic ACDM in 
the high-redshift regime as well as produce ACDM ac- 
celeration of the universe at low redshift without a true 
cosmological constant; it should also fit Solar system ob- 
servations. 

The particular form chosen by [23] is: 



Ci{R/'m,^ 



C2(i?/m2)" + 1 



with 



m 



(8315Mpc)" 



0.13 



(2) 



(3) 



where k = 8nG and po is the average density today. In 

this model, modifications to GR only appear at low red- 
shift, when we are safely in the matter dominated regime. 
The properties of the model are well described by the 
auxiliary scalar field ff, 



dJ{R) 
dR 



Before going to the expansion history, it is worth briefly 
reviewing the main features of this model, following the 
original presentation in [23]. The trace of the modified 
Einstein equations serves as the equation of motion for 
fR-- 



3n/fl-i? + /fli?-2/^ 



(4) 



and thus is essentially the Compton wavelength of fa at 
the background curvature in units of the horizon length. 

In the static limit with <C 1 and |//i?| <C 1, Eqn. 
4 becomes 



(9) 



where p is the local density. This equation has 2 modes 
of solutions. One is the very high curvature R k, n^p 
and the other one is the low curvature (but still high 
compared to the background density) R <C K^p. For more 
on the interplay of these two regimes and applications in 
solar system observations sec [23]. 

Let's now move on to the expansion history. For the 
model to yield behavior that is observationally viable re- 
quires a choice for the present day value of the field 
fm ^ 1- This is equivalent to Rq^ m^. In that case the 
approximation R^ is valid for the whole expansion 
history and we have: 



lim f{R) « m + -^m —— 

m^/R^O C2 (4 \ R 



(10) 



In the limiting case of ci/c^ ^ at fixed ci/c2 we ob- 
tain a cosmological constant behavior ACDM. Thus to 
approximate the ACDM expansion history with a cos- 
mological constant Qa and matter density flm we set: 



Cl 
C2 



6 



_0a 



(11) 



This leaves 2 remaining parameters to play with: n and 
C1/C2 to control how closely the model mimics ACDM. 
Larger n mimics until later in the expansion history, while 
smaller ci/c^ mimics it more closely. 

For flat ACDM we have the following relations: 



or, in terms of the effective potential, 



(5) 



The effective mass for the Jr field is then given by the 
second derivative of Vcs, evaluated at its extremum: 



(12) 



(13) 



2 1 fi + fR p 

mgff = -^TTo- = - I — R 



dfl 3 



Irr 



(6) 



The Compton wavelength of the field is then given by 

^/« = "^^ff • 

It is very convenient to introduce a dimensionless quan- 
tity: 



p _ fRR n' H 



(7) 



It has been shown [23] that in the high-curvature regime 
B is connected to the Compton wavelength via: 



(8) 



As we will see later these are the necessary ingredients 
for the application of perturbation theory to the model. 

Finally we need to obtain a suitable parametrization 
of the mode. At the present epoch we have: 



Ro 



Cl 



RO 



12 

12 
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(14) 



(15) 



In particular, for 
Ro = Aim? and Jrq « 



0.76 and Q„i = 0.24, we have 
-n Ci/c2/(41)"+^. ^From now 
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on we will parametrize the model through f^o n. 
Fig. 9 in [23] shows that there is a wide range of viable 
parameter values which satisfy Solar system and Galaxy 
requirements. We re-iterate that what makes this par- 
ticular f{R) model viable is its Chameleon behavior - 
the possibility of uniting galaxy and solar system obser- 
vations with the expansion of the Universe. 



III. PERTURBATION FORMALISM 

By definition, the dark sector (dark matter and dark 
energy) can only be inferred from its gravitational con- 
sequences. In general relativity, gravity is determined by 
the total stress-energy tensor of all matter and energy 

We may consider the Hubble parameter H{z) to be 
fixed by observations. In a dark energy model, p is given 
by the Friedman equation of GR: p = SH^/SttG. The 
equation of state parameter is w = — 1 — 2H/3H'^. The 
corresponding modified gravity model has matter density 
to be determined from its Friedman-like equation. We 
will consider MG models dominated by dark matter and 
baryons at late times. 



A. Metric and fluid perturbations 

With the smooth variables fixed, we will consider per- 
turbations as a way of testing the models. In the Newto- 
nian gauge, scalar perturbations to the metric are fully 
specified by two scalar potentials and 



ds'^ = -(1 - 2*) dt^ + {1- 2$) a2(t) df^ 



(16) 



where a{t) is the expansion scale factor. This form for 
the perturbed metric is fully general for any metric the- 
ory of gravity, aside from having excluded vector and 
tensor perturbations (sec [24] and references therein for 
justifications). Note that corresponds to the Newto- 
nian potential for the acceleration of particles, and that 
in general relativity $ = — ^ in the absence of anisotropic 
stresses. 

A metric theory of gravity relates the two potentials 
above to the perturbed energy-momentum tensor. We in- 
troduce variables to characterize the density and velocity 
perturbations for a fluid, which we will use to describe 
matter and dark energy. The density fluctuation d is 
given by 



S{x,t) 



p{x, t) - p{t) 

m 



(17) 



where p(x, t) is the density and p{i) is the cosmic mean 
density. The second fluid variable is the divergence of the 
peculiar velocity 



where v is the (proper) peculiar velocity. Choosing 
instead of the vector v implies that we have assumed v 
to be irrotational. Our notation and formalism follows 
that of [25]. 

In principle, observations of large-scale structure can 
directly measure the four perturbed variables introduced 
above: the two scalar potentials \E' and and the density 
and velocity perturbations specified by 5 and 6. It is 
convenient to work with the Fourier transforms, such as: 



/ 



5{k,t) = / d^x S{x,t) e 



-ik-x 



(19) 



When we refer to length scale A, it corresponds to a 
a statistic such as the power spectrum on wavenumber 

k = 2tt/\. Wc will henceforth work exclusively with the 
Fourier space quantities and drop the " symbol for conve- 
nience. 



B. Linearized fluid equations 

We will use the perturbation theory equations for the 
quasi-static regime of the growth of perturbations. Wc 
begin with the fluid equations in the Newtonian gauge, 
following the formalism and notation of [26] . 

For minimally coupled gravity models with baryons 
and cold dark matter, but without dark energy, we can 
neglect pressure and anisotropic stress terms in the evo- 
lution equations to get the continuity equation: 



- - 3$ ~ — 
a I a 



(20) 



where the second equality follows from the quasi-static 
approximation as for GR. The Euler equation is: 



-H9- 



(21) 



We parametrize modifications in gravity by two func- 
tions Gcff{k,t) and ri{k,t) to get the analog of the Pois- 
son equation and a second equation connecting $ and 
^ [25, 27]. We first write the generalization of the Pois- 
son equation in terms of an effective gravitational con- 
stant GeS- 



A;2$ = -4'7rGes{k,t)pMGa^SMG 



(22) 



Note that the potential $ in the Poisson equation comes 
from the spatial part of the metric, whereas it is the 
"Newtonian" potential ^ that appears in the Euler equa- 
tion (it is called the Newtonian potential as its gradient 
gives the acceleration of material particles). Thus in MG, 
one cannot directly use the Poisson equation to eliminate 
the potential in the Euler equation. A more useful ver- 
sion of the Poisson equation would relate the sum of the 
potentials, which determine lensing, with the mass den- 
sity. We therefore introduce GeS and write the constraint 
equations for MG as 



p) = V-v, 



fc2($-*) 



-87rGeff(A:, t)pMGa^SMG 



(23) 
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FIG. 1: Fractional change at z = in the density power spectrum 
of the f{R) model compared to AC DM for a set of choices of fno- 
The upper panel is the n = 4 model, while the lower panel has 
n = 1. This figure can be compared with Fig. 4 in [23]. 



FIG. 2: Fractional change at 2 = in the velocity power spectrum 
of the f{R) model compared to ACDM for a set of choices of fno- 
The upper panel is the n = 4 model, while the lower panel has 
n = 1. 



$ = -^1 ri{k, t) 



(24) 



where GcS = Gcff (1 + ??"^)/2. 

The parameter Gcff characterizes deviations in the ($— 
"^yS relation from that in GR. Since the combination 
<i> — 5" is directly responsible for gravitational lensing, 
Gcff has a specific physical meaning: it determines the 
power of matter inhomogeneities to distort light. This 
is the reason we prefer it over working with more direct 
generalization of Newton's constant, Gcff. 

With the linearized equations above, the evolution of 
either the density or velocity perturbations can be de- 
scribed by a single second order differential equation. 
From Eqns. 20 and 21 we get, for the linear solution, 

initial 

{k)D{k,t), 



S + 2H5 



= 0. 



(25) 



For a given theory, Eqns. 23 and 24 then allow us to 
substitute for in terms of 5 to determine D{k,t), the 
linear growth factor for the density: 



D + 2HD 



pa^ D = 0. 



(26) 



What we need now are Gcfi and rj for our modified 
gravity model. 



C. Parametrized post-Friedmann framework 

Hu & Sawicki [23] also developed a formalism for si- 
multaneous treatment of the super-horizon and quasi- 
static regime for modified gravity theories (in particular 
f{R) and DGP) [28]. They begin by describing the differ- 
ent regimes individually and the requirements they im- 
pose on the structure of such models. Consequently they 
describe a linear theory parametrization of the super- 
horizon and the quasi-static regime and test it against 
explicit calculation. What is important for this paper is 
the proposed interpolation function for the metric ratio 



9 = 



$ + ^ 
$ - ^' 



(27) 



For consistency we will express formulae from Hu & 
Sawicki in terms of physical time instead of expansion 
factor. We start with a background FRW universe for 
which we have the curvature in terms of the Hubble 
parameter R = 6H + 12H^. As we mentioned the 
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background evolution H is chosen to match that of flat 
ACDM. We can then compute the Compton parameter 
B for our preferred model since we know what f{R), R 
and H are: 



The next step is to look at the super-horizon regime. 
In that case we know how to calculate the potentials $ 
and ^. Their evolution is given by [13]: 



«>+ 1 



H 



B 



HH H{l-B) m 
H H B 



B^ 1 



'H^ HH H{l-B) 



H^^ = 0, (kH 0) 



-^-B^/H 

* = — ,(fcg^o) 



(29) 



(30) 



This allows us to compute gsH = g(t, ku — > 0). 

Furthermore in the case of subhorizon evolution where 
ku = k/a{t)H » 1 we have gqs = —\ [13]. 

According to [28] in the case of f{R) theories we can 
use the following interpolation functions for the metric 
ratio: 

9SH + gQs{.CgkHT^ 



g{t,k) = 



1 + [cgknY^ 



(31) 
0.71^1/2 



The evolution then is well described by [28] Cg 
and rig = 2 In terms of the post-Newtonian parameter: 
r] = —^/^ wc have g = {rj — l)/(r7 + 1). 

We still need one more ingredient, Geff , which is given 
by [28] 

Geff(fc,i) = -^ (32) 

1 + In 

We now have all the needed components to use in equa- 
tion (26) for the growth factor of density perturbations. 
Calculations of the fractional change in the linear density 
power spectrum compared to GR have been done [23]. 
We show these in Fig. 1, using Geff and r] for the partic- 
ular model wc investigate. We have assumed the transfer 
function for the concordance A— CDM model consistent 
with the 5-year WMAP data [29]. 

We can also use the relations given above to obtain the 
linear growth factors for 6 and the potentials from D. 
Note that in general the growth factors for the potentials 
have a diS'erent k dependence than D. 



IV. THE BISPECTRUM IN PERTURBATION 
THEORY 

A. Nonlinear Fluid Equations 

The fluid equations in the Newtonian regime are given 
by the continuity, Euler and Poisson equations. Keeping 



the nonlinear terms that have been discarded in the study 
of linear perturbations above, the continuity equations 
gives: 



6 + 



-I 



d^ki k- ki 
(27r)3 fcf 



0{ki)S{k - ki) (33) 



where the term on the right shows the nonlinear coupling 
of modes. Note that the time derivatives are with respect 
to conformal time in this section. 
The Euler equation is 



e + H0 + 



, 2 T 

J (27r. 



)^ 2fc?|fci -fcsl^ 



(34) 

We neglect pressure and anisotropic stress as the energy 
density is taken to be dominated by non-relativistic mat- 
ter [30]. The Poisson equation is given by Eqn. 23 and 
supplemented by the relation between \1/ and $ given by 
Eqn. 24. Using these equations we can substitute for 
in the Euler equation to get 



9 + He 



SttG 



off 



(1 + ^7) 



pMca 5 



(27r)'^ 2kf\ki - fc2p 



Eqns. 33 and 35 are two equations for the two variables 
5 and 0. They constitute a fully nonlinear description 
and can be solved once r) and Gcff are specified. An 
important caveat is that they may nevertheless be in- 
valid on strongly nonlinear scales or for particular MG 
theories. For the f{R) model considered here, they are 
valid on scales well above 1 Mpc; on smaller scales the 
chameleon mechanism modifies the growth of structure 
[31]. Since we will use perturbation theory, our approach 
breaks down once ^ ~ 1 in any case. 

Next we consider perturbative expansions for the den- 
sity field and the resulting behavior of the power spec- 
trum and bispectrum. Let 6 = -|- ^2 + ••• where 
62 ~ 0{Si). In the quasilinear regime, i.e. on length 
scales between ~ 10-100 Mpc, mode coupling effects can 
be calculated using perturbation theory. While this is 
strictly true only for general relativity, a MG theory that 
is close enough to GR to fit observations can also be ex- 
pected to have this feature. 

For MG, following [25], let us simplify the notation by 
introducing the function: 



ak,t) 



(36) 



which is simply 47rG in GR but can vary with time and 
scale in MG theories. The evolution of the linear growth 
factor is given by substituting for "if in Eqn. 37 to get 
(as above, but with conformal time here) 



61 + H61 - Cpa^Si = 0. 



(37) 



6 



The linear growth factor has a scale dependence for our 
f{R) as shown in Fig. 1. 

In addition we show below that the second order so- 
lution has a functional dependence on Gcff and rj that 
can differ from GR. Thus potentially distinct signatures 
of the scale and time dependence of Gcs{k,z) can be 
inferred from higher order terms. These rely either on 
features in k and t in measurements of and Ps-, or 

on the three-point functions which can have distinct sig- 
natures of MG even at a single redshift [32] . Quasilinear 
signatures due to r]{k, z) can also be detected via second 
order terms in the redshift distortion relations for the 
power spectrum and bispectrum. Our discussion general- 
izes that of [6] who examined a Yukawa-like modification 
of the Newtonian potential. 



B. Second order solution 

From a perturbative treatment of Eqns. 33 and 35 the 
second order term for the growth of the density field is 
given by [25] 



Hh[5x,5{\+l2[5iJi] 



A[<5'i,<5i], 



(38) 



where Ii and I2 denote convolution like integrals of the 
two arguments shown, given by the right-hand side of 
equations 33 and 35 as follows 



/i[<5i,<5i](fc) 



^t^6,{k,)5,{k-k,) (39) 



and 



h[Su5i]{k) 



(27r)3 kl 
(fiki k% • (fc- fci) 

(2^)3 2fc2|fci-fc2|2 



6i{ki)Si{k ~ kl). 

(40) 

Finally, the last term in Eqn. 38 is simply Ii[6i,6i] = 
-^i[^i,<^i] + -^i[<^i, i^i]- Note that by continuing the itera- 
tion, higher order solutions can be obtained. 



C. Three-point correlations 

Distinct quasilinear effects are found in three-point cor- 
relations - we will use the Fourier space bispectrum. Re- 
cently Tatekawa & Tsujikawa have performed a similar 
perturbative analysis and presented results on the skew- 
ness [34]. We prefer to use the bispectrum as it allows 
one to study specific quasilinear signatures contained in 
the configuration dependence of triangle shapes. The bis- 
pectrum for the density field Bs is defined by 

{S{ki)S{k2)S{k3)) = {2TrfSuiki + k2 + k3)Bs{h,k2,k3) 

(41) 

Since Bg ~ (6^) ^ (5f52) (using {6f) =0 for an initially 
Gaussian density field) at tree level, the second order 
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FIG. 3: Upper panel: The Bispectrum of the f{R) model for 
equilateral triangles depending on scale for the f{R) model with 
/-RO = 10~^, n = 1. The corresponding regular gravity bispectrum 
is shown as the dashed curve. Lower panel: The Ratio of the f(R) 
Bispectrum to the regular gravity one. 



solution enters at leading order in the bispectrum. Note 
also that the wavevector arguments of the bispectrum 
form a triangle due to the Dirac delta function on the 
right-hand side above. 

The bispectrum is the lowest order probe of gravita- 
tionally induced non-Gaussianity. A useful version of it 
is the reduced bispectrum Q, which for the density field 
6 is given by 



Bs{ki,k2, fcs) 



Ps {ki)Ps (fcs) + Ps {k2)Ps ih) + Ps {kl )Ps (fca) ' 

(42) 

Q is useful because it is insensitive to the amplitude of 
the power spectrum; thus e.g. at tree level and in the 
case of a scale free linear power spectrum P{k) ~ /c" it 
is static and scale independent for regular gravity[33]. 



V. RESULTS 

In Fig. 3 we present the bispectrum for the f{R) model 
and its dependence on scale for two different redshifts. 
For comparison we also show the bispectra predicted in 
GR. In the lower panel we present the ratio of the f{R) 
bispectra to that in GR. We have assumed the transfer 
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FIG. 4: The reduced bispectrum Q for equilateral triangles for the 
f{R) model compared to regular gravity. Top panel: fno = 10~^, 
n = 1. Middle panel: fno = 10"'', n = 4. Bottom panel: fug = 
10-5, n = 4. 



function for the concordance A— CDM model consistent 
with the 5-year WMAP data [29] . Calculation is done at 
tree level with fi™ = 0.24. The bispectra in f{R) gravity 
are enhanced relative to GR, increasingly so at high-fc. 
The enhancement is of order 10-20% for observationally 
relevant scales around fc ~ 0.1 and redshifts below unity. 

We turn our attention to the reduced bispectrum Q, 
which is expected to show features not associated with 
the linear power spectrum [33]. We show two relevant 
cases. The first one is with equilateral triangles, shown 
in Fig. 4. Deviations from GR are at the percent-level, 
which makes them impossible to detect with current mea- 
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FIG. 5: The reduced bispectrum Q for isosceles triangles (ratio 
of sides lengths 1:3:3), where the x-axis shows the length of the 
smallest k in the triangle. The ratio of Q for the f{R) model with 
Jro = IQ— 5, n = 1 to regular gravity is shown. 



surements. Qualitatively, the parameters of the model do 
influence the scale and time dependence of the reduced 
bispectrum. 

Slightly stronger deviations from regular gravity are 
observed for isosceles triangle configurations, shown in 
Fig. 5. Once again strong scale and time variation is 
observed when changing the model parameters. We also 
show the angular dependence of the reduced bispectrum 
for 3:1 ratio configurations and its comparison to regular 
gravity in Fig. 6. 



A. Does the Linear Growth Factor Determine 
Nonhnear Clustering? 

Quasilinear effects, and the bispectrum in particular, 
show signatures of gravitational clustering. However, we 
find that the reduced bispectrum Q remains very close 
to to that of GR in the modified gravity models we have 
considered (see also [34] for a related study). This is 
qualitatively similar to previous findings about the in- 
sensitivity of Q to cosmological parameters within a GR 
context [33]. 

Thus the deviations in the bispectrum are largely de- 
termined by the linear growth factor. Ambitious future 
surveys will be needed to achieve the percent level mea- 
surements required to probe the unique signatures of 
modified gravity in the bispectrum. On the other hand, 
this means that the bispectrum can be used as a consis- 
tency check on the power spectrum. It has been shown 
that the bispectrum contains information comparable to 
the power spectrum, thus improving the signal-to-noise 
[35] . Equally importantly, it is affected by sources of sys- 
tematic error in different ways and contains signatures of 
gravitational clustering that are unlikely to be mimicked 
by other effects. 

Currently the results of N-body simulations for this 
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FIG. 6: Upper panel: The angular dependence of the reduced 
bispectrum Q for triangles with the two fixed side lengths in the 
ratio 1:3 and k = 0.1/i/Mpc for the smaller one. The f{R) model 
with fnQ = 10~^, n = 1 at 2 = is used. The dashed line shows 
the prediction for GR. Lower panel: Comparison to regular gravity 
at two redshifts. We can observe that for obtuse shapes (6 ~ 0) our 
model predicts a lower Q for f{R) gravity, while for shapes close 
to isosceles it is enhanced. 



f{R) studies of [31] and [40] show distinct effects of the 
chameleon field on small scales (comparable to galaxy 
and cluster sized halos for most models); [28] suggest a 
fit for the nonlinear power spectrum with additional pa- 
rameters that describe the transition to the small-scale 
regime. The DGP study of [38] requires inclusion of non- 
linear terms in the Poisson equation. So clearly for differ- 
ent models there can be new nonlinearities and couplings 
to additional fields that impact the small-scale regime of 
clustering. Even so, for a class of models that includes 
the f{R) models studied here, the quasistatic, Newtonian 
description parameterized by Geff and rj = ^/'i! applies 
over a wide range of scales relevant for large-scale struc- 
ture observations. In this regime, to a good approxima- 
tion, many clustering statistics can be predicted using the 
linear growth factor and the standard gravity formalism. 



B. Implications for Lensing and Dynamics 




0.005 0.010 0.050 0.100 0.500 
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f(R) model [31] show that our assumptions are consistent 
with their result for the power spectrum in the regime of 
validity of perturbation theory (e.g. for k < Q.3[h/Mpc]). 
On smaller scales, the chameleon mechanism modifies the 
power spectrum. 

Modifications to the Newtonian potential were simu- 
lated by [20] (also see [6], [36]). These simulation stud- 
ies found that, to a good approximation, the nonlinear 
power spectra depend only on the initial conditions and 
linear growth. The standard fitting formulae for Newto- 
nian gravity [37] were adapted to predict the nonlinear 
spectrum at a given redshift. Therefore tests for mod- 
ified gravity using the power spectrum require either a 
measurement of the scale dependent growth factor (in 
combination with the initial spectrum measured from the 
CMB), or of measurements at multiple epochs. More 
likely a combination of probes will be used for robust 
tests of gravity (see e.g. [25]). 

Related studies of nonlinear clustering in f{R) models 
or DGP gravity are in progress [38-40] ; these authors are 
considering the power spectra, bispectra as well as halo 
properties, in particular the halo mass function. The 



FIG. 7: The fractional change in r (blue lines) and /3-^ (red lines) 
for the f{R) model with /^g = 10~^, n = 1 compared to GR as a 
function of scale. 

Lensing observations provide estimates of the conver- 
gence power spectrum and bispectrum (see e.g. [41]). For 
a rough estimate of these quantities, we take the source 
galaxies distribution to be a delta function at a given red- 
shift and take the lensing matter to be situated at half 
the distance. Since we have taken the expansion history 
to be ACDAI, so comoving distances are the same as in 
GR. With these approximations, From Eqn. 

23 we can see that the difference of the lensing behavior 
of our modified gravity case and the regular gravity one 
is given by (see also [32]) 

p p / Goff(Ml/2) V 

P^-^s, cx Ps I — I (43) 

where ti/2 is the physical time at the redshift of the lens- 
ing matter. 
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Thus we can write: 



PuMG PSMG 
(X 

PkGr Psgr 



G 



(44) 



Analogously we can see that the ratio of convergence 

bispcctra behaves like {Gca{k,ti/2)/G)'^ but the ratio of 
the reduced bispectra behaves like 



QkMG Qsmg 

(X 



G 



Qk.gr Qsgr Gcs{k,ti/2) 



(45) 



Thus the convergence power spectrum and reduced 
bispectrum could be successfully used to differentiate 

and/or rule out models of modified gravity. Unfortu- 
nately in the currently discussed model we have Eqn. 32; 

r(fc,^(a)) = ^= ) (46) 

which deviates from unity at the order of fm, which is 
much smaller than unity. Thus the convergence power 
spectrum and reduced bispectrum follow almost identi- 
cally the predictions for their matter counterparts. This 
means that the comparison of lensing to tracers of mass 
fluctuations does not reveal distinct signatures of f{R) 
gravity. 

The Newtonian potential ^ drives dynamical observ- 
ables such as the redshift space power spectrum of galax- 
ies [25], [27]. The velocity growth factor is related to 
the density growth factor via the function /?: 



De oc apHD 



(47) 



where f3 = d\n D/dlna. This function varies with scale 
and expansion factor for our f{R) models. Both t and /? 
can be seen in Fig. 7 which clearly shows that observables 
based on peculiar velocities would show a clear signature 
of f{R) gravity. More detailed calculations of lensing and 
velocity statistics are beyond the scope of this study. 



three-point correlation fimction, the bispectrum. of mat- 
ter. Our results are applicable up to scales at which the 
derivation of the quasi-linear perturbation theory is still 
valid. 

We found that in the bispectrum the dominant behav- 
ior is due to the difference in the linear growth factors 
between the modified gravity model and regular grav- 
ity. The bispectrum itself shows significant departures 
in scale and time compared to the predictions of GR. 
However the reduced bispectrum, which is independent 
of the linear growth factor in perturbation theory for GR, 
remains within a few percent of the regular gravity pre- 
diction. It does show interesting signatures of modified 
gravity at the percent level. 

Our results are consistent with studies of the nonlinear 
regime via N-body simulations, which have found that on 
a wide range of scales the nonlinear power spectrum can 
be predicted using the (modified gravity) linear growth 
factor in the standard formulae developed for Newto- 
nian gravity. Our results imply that three-point corre- 
lations follow this trend at the few percent level. (The 
regime around and inside halos probed by [31] to test the 
chameleon behavior is not included in our perturbative 
study.) It would be interesting to compare perturbative 
and simulation results for the bispectrum for the models 
considered here and other modified gravity models. 

Upcoming surveys in the next five years will not attain 
the percent level accuracy at which the reduced bispec- 
trum shows distinct signatures of f{R) gravity. In this 
time- frame, the bispectrum will be useful as a consistency 
check on potential deviations from GR found in the power 
spectrum. Such a check is useful since measurement er- 
rors and scale dependent biases of tracers can mimic some 
of the deviations in the power spectrum. Next genera- 
tion surveys, to be carried out in the coming decade, will 
provide sufficient accuracy to test the distinct signatures 
seen in our results for the reduced bispectrum. With 
these surveys, the bispectrum can provide truly new sig- 
natures of gravity. 



VI. DISCUSSION 

In this paper we studied the growth of structure in 
an f{R) modified gravity model [23]. In the quasilin- 
ear regime, we used perturbation theory to calculate the 
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